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Abstract 

The Ito and Stratonovich approaches are carried over to quantum stochastic systems. Here the white noise rep- 
resentation is shown to be the most appropriate as here the two approaches appear as Wick and Weyl orderings, 
respectively. This introduces for the first time the Stratonovich form for SDEs driven by Poisson processes or quan- 
tum SDEs including the conservation process. The relation of the nonlinear Heisenberg ODES to asymptotic quantum 
SDEs is established extending previous work on linear (Schrodinger) equations. This is shown to generalize the clas- 
sical integral transformations between the various forms of stochastic calculi and to extend the Khasminskii theorem 
to the quantum setting. 
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r— nl Introduction 

stochastic asymptotic analysis of dynamical systems has its origins in Einstein's theory of Brownian motion, however, 
i-£h it extends to very general classes of systems, both classical and quantum. 

in the classical case, one considers typically a dynamical variable x t = Xt{\) determined from an ODE -^xt{\) = 
\F(xt{\), Mt), xq(X) = xo, where A > is a coupling parameter, M t is an external (stochastic) input and F is some 
J function (usually Lipschitz with bounded first partial derivatives). Supposing that ^ J^ +t ds F(x, M s ) converges as 
T — > oo in probability to F {x) uniformly in t for arbitrary x, then under suitable conditions [1] the averaging principle 
^ states that x t /\(\) converges in probability to Xt, the solution of the averaged ODE -^Xt = F (x t )', %o = %o- One may 
^hink of x t /\(X) as a perturbation away from the averaged solution. If F Q = 0, the averaging principle states that x t /\(\) 
will not vary significantly from xq on time scales [0, T/X] [2]. This lead Stratonovich in 1961 to suggest that on time scales 
[0, T/A ] stochastic fluctuations accumulate appreciably; this idea was substantiated by Khasminskii [3] who proved that 
£NJ yj/A 2 (A) converges weakly to a Markov diffusion. In this context, the interpretation in terms of a white noise Langevin 
equation, originally given by Wong and Zakai [4], has been revealing. If £(i) are regular stochastic processes with mean 
zero and correlation E[£(i)£(s)] = yj-c(^-), so that the A — >• limit leads to a white noise process, then the dynamical 
equations of motion can converge to an asymptotic Langevin equation which retains the pre-limit form provided one uses 
• — the Stratonovich version of stochastic calculus [5] . 

When dealing with quantum systems the scaling employed in the Khasminskii theorem corresponds to the van Hove 
^ re-scaling of the dynamical variables [6] . The formulations of most quantum Langevin equations encountered in physics 
"[7] can occasionally be justified as an asymptotic quantum stochastic limit of such re-scaled variables [S]. The germ 
of the approach to be developed in this paper comes from the work of Accardi, Frigerio and Lu [9] wherein a class 
of noncommutative Khasminskii theorems are established. The point of view presented here is closer in spirit to that 
of Wong and Zakai. In particular, quantum white noises [10] play a fundamental role. These are Bose operators (t) 
satisfying a canonical commutation relation (CCR) 

[a(t),a{s))] =K8+(t-s) + K*S-(t-s), (1.1) 

where k is a complex number and S±(t) are special delta functions on the space of regulated functions picking out the 
future/past values at a time t. The form (1.1) arises as the limit form of the CCR for physical fields wherein the right- 
hand side would be typically the causal (Feynman) propagator. The a"(i) have white noise spectra and play a role similar 
to the input processes studied by Gardiner [11], though with different operational properties and without the effective 
restriction k = | . As such, the choice between Ito and Stratonovich versions of stochastic calculus corresponds to the 
choice between Wick and Weyl ordering, respectively, of the noises with respect to the integrands. Note that the product 



1 



of two Wick ordered expressions will require further re-ordering to attain Wick order again, this leads to extra terms 
which correspond to the standard shift accounted for by the quantum Ito table [12]. This notion, in fact, dates back to an 
observation of Hudson and eater [ 131. It should also be mentioned that an alternate description of Stratonovich quantum 
stochastic calculus, not involving white noise and taking k = h, has been given independently by Chcbotarev [14]. It is 
worth mentioning that while the stochastic limits described here should preserve canonical structures, especially unitarity, 
many important physical properties such as detailed balance and the KMS condition are typically lost [15]. Heuristic 
accounts of these white noise processes have appeared [16, 171 along with an account of their applications to physic [18]. 
However, only linear systems have been treated up to now. In this paper, the analysis is extended to nonlinear systems 
which essentially emerge as Heisenberg evolutes corresponding to linear quantum SDEs (i.e. stochastic Schrodinger 
equations). Complete transformation laws are derived. The classical formulae for Ito-Stratonovich conversion for both 
Wiener and Poisson driven systems are recovered as special cases. 



2 Classical stochastic differential equations 



In this section several well-known features [19] of SDEs are reviewed. Let it, er^ <G C(™xM, W 1 ) satisfy appropriate 

Lipschitz and growth conditions, and let m[°^ be a set of stochastic processes on a probability space (CI, P). Let 
x e L 2 (fl, T 

equation (2.1) below is technically meaningless, 



be independent of the processes m\°^ for t > 0. Whenever the are nondifferentiable, the formal 
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X = Xq. 



(2.1) 



Different finite-step numerical schemes used to define an estimate for X, in (2.1) lead to inconsistencies in the small 
time-step limit which do not arise when dealing with ODES. For instance, let V be an ordered partition of [0,i], say 
[0 = to < t\ < ■ ■ ■ < t?j = t], and denote max|ij +1 — t 3 \ by \V\. One constructs an approximate solution Xl^ to (2.1) in 
the form of a random variable Xf N , obtained by the iterative scheme 
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One then understands the limit of X^ for finer partitions as a mean-square limit, 
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J* |« (X s ,s) ds + £ ct(q) (*a, s) dAfW | := ^XT N 



(2.2) 



(2.3) 



if such a limit Xf G L 2 (fi, J 7 , P) exists. Amongst the various schemes to compute 

(2.4) 



the increment AX^, the simplest is the Eulcr method 

AX% := v (X^tj) [t 3+1 t 3 ] + £ fe, t 
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an alternative is given by the averaging scheme 

: =\[v (X? j+1 , t j+1 ) + v (Xf., tj )] [t j+1 - tj 
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(2.5) 



If the were zero, so that (2.1) reduces to an ODE, then the Euler scheme has local error of 0(h 2 ) : whereas the 
averaging scheme has local error of 0(h 3 ), where h is the largest increment size (which is, of course, proportional to IPI 
for the deterministic case). If the a'"' are nonzero then the noise increments M tj+1 — M tj may fluctuate rapidly enough 
to lead to h = 0(\V\ a ) with a < 1. For instance, in the case of the Wiener process as noise one has a = h. In this 
situation the Euler scheme leads to an ( I? 13 1 ) local error which makes a nontrivial contribution to the stochastic integral 
(2.3) which is not present in the averaging scheme. The respective solutions (when they exist) are called the Ito and 



Stratonovich stochastic integrals and the following notation shall be used: X t = Xq + f* dX s , with 

(Ito) dX t = v(X t ,t)dt + J2<T {a) (X u t)dM$ a \ (2.6a) 

a 

(Stratonovich) dX t = v {X u t)dt + ^ cK Q ) (X t ,t)dM^ a) . (2.6b) 

a 

The deterministic component of the integrand (coefficient of dt) is taken to be Ricmann integrable and so the over- or 
under-bar notation can be omitted, as frequently will be the case. In general, for X t and Y t stochastic processes, one 
understands the limits X t ., = 0. 



,i N-l t _ N-l 

/ X s dY s := l.i.m. V XT AY f v , / X s dY s := l.i.m. V XfAY t p , (2.7) 
Jo VHof^ <> *»' J l*H>°^ 3 3 



where AY P := Y^ +i — Y p and Xf. := | \XT j+1 — -^^J ■ Moreover, to compute the increments of the product X, Y, one 
notes 

A (XY) V = X v Y v — X V Y V = X V AY V + AX V Y V + AX V AY V 

= ~Xj.AYf. + AX?Y£. (2.8) 

The terms AX AY need not be negligible; for the Wiener process it is 0(|"P|). Thus, if the noise processes M t (Q) admit 
a stochastic calculus and if X t ,Y t are processes driven by these noises, then one expects that the Ito form breaks the 
Leibniz rule of calculus, 

d (X t Y t ) = X t dY t + dX t Y t + dXtdY t , (2.9a) 

while the Stratonovich form retains it 

d {X t Y t ) = XtdYt + dX t Y t ; (2.9b) 

here the equivalence relation means equality up to O(dt) in the deterministic part which implies equality under the 
integral sign. Using (2.9b) inductively, one sees that 

dX 2 t = 2X t dX tl dX? = {ixf + Xf) dX t , dX? = (xf + %xf + 2xf) dX t , ■■■ . (2.10) 
Now Xj = | (X t + dX t ) n + |X t ", from which it follows that Xf - X~t = O UdXtf} and so 

XfdXt = xTdX t + o ((dX t ) 2 ) • 

In the case of diffusion processes, one then has X™dX t = X t n dX t , and so dX™ = nl t ™ 1 dX t which implies that the 
chain rule formula of standard calculus holds when applying the Stratonovich calculus to ditfusions. However, one can 
see already that this will not be true for processes driven by the Poisson process with nonlinear noise coefhciciits. 

Examples Whereas the Stratonovich integral with respect to general martingales can be formulated, it is interesting to 
give explicitly the Wiener and Poisson cases. 

Wiener Process as noise: The Wiener process B tl admits a stochastic calculus with (dB t ) 2 = dt and higher powers 
o(dt). The Stratonovich SDE dX t = v(X t ,t)dt + <j(X t ,t)dB t is well known to be equivalent to an Ito SDE dX t = 
v(Xt, t)dt + a(Xt, t)dBt, where the coefficients are related by 

V = V + 2°' (T '' cr = £7, (2-11) 

with a'(x,t) := J^cr(x,t). Inversely, v = v — \oa. 

Poisson process as noise: The Poisson process Nt admits a stochastic calculus with (dNt) n ~ dN, for all integers 
n >i. The fluctuations in dNt are 0(dt). The Stratonovich SDE dXt = v(Xt,t)dt + /j,(Xt,t)dNt, is equivalent to the Ito 
SDE dX t = v(X t ,t)dt + fl(X t ,t)dN t , where the coefficients are related by 

v = v, p,(x,t) = - [ii(x + {t(x,t)) +fj,(x,t)] . (2.12) 



The Stratonovich form for the Poisson integrals is next derived explicitly. Let / = fix) be analytic, then by Taylor 
expansion one has 




Therefore, multiplying by AiV tj = N tj+1 — N tj , and re-summing leads to 

/ AN tj = f (xf. + p, (X*, f 3 -)) AN t] 

On account of this, 

= ^+l{v(xr j+1 , tj+l )+v(xr j ,t j )}At j 

+ \{^ (X% +1 , t 3+1 ) + n (Xf., tj ) } AN tj 
= Xf. + v (Xf.,t 3 ) At, + 1 { M (Xf. + jj, (X%, tj) , tj) + M (X%, tj) } AN tj . 
and so the relation (2.12) closes the transformation. 

Remarks i) In the case of linear Poisson-noisc coefficient fi(x) = fi x + c the relation (2.12) gives jl{x) = (n x + c)/(l — 
i/i ). The requirement on the coefficients is that x >— > x + n(x, t) is strictly monotone, for each t, and so defines a proper 
change of variable. For general ^(x) there exists no closed expression for /2, however, it may be approximated by the 
iteration p, n+1 (x) = |{/i(x + fi n {x)) + /J,{x)}, fi (x) = fi(x). 

ii) There are other numerical schemes which can be employed. For instance, the mid-point 2nd order Runge-Kutta 
method 

xf ]+i =xT ]+ v (.y/\,:) Atj + £ ^ (-vr./;) am<°>, 

where t* := ^ (tj+i + tj) and xjl is obtained from an Euler time-step from tj to £* 

xT* = x% + v (a/:./:) [t* - tj] + £ - {a) 1*4? ~ < J ] 

a 

with v, the associated Ito coefficients. This procedure is, in fact, equivalent to the Stratonovich choice. However, the 
method given by is equivalent to the Stratonovich choice only for the Wiener process as noise. For the Poisson process, 
it coincides with the Stratonovich choice only if /j,(x, t) is linear in x; otherwise it leads to a stochastic integral equivalent 
to an Ito integral with coefficient fi(x,t) related by n{x + ^fi(x,t),t) = fi(x,t). 

iii) In general X t is said to be Riemann integrable with respect to Y t , if the mean-square limit of JA Xf. AY^ exists 
whenever an arbitrary prescription to determine s G [tj,tj + i] is given. The Ito calculus is based on the choice Sj = tj. 
One may call this the retarded Ito theory. The advanced Ito theory can be defined via the prescription that Sj = 
tj+i- The Stratonovich calculus then corresponds to the averaging over the retarded and advanced Ito versions. An 
alternative definition of Stratonovich integrals as simply mean-square Riemann integrals with Sj = t* is frequently taken 
in probabilistic literature, however, as seen in the previous remark this does not always coincide with the definition here. 

3 Quantum Stochastic Calculus 

We begin by recalling the Hudson-Parthasarathy theory of noncommutativc stochastic processes [12, 20]. Let T-Lq and 
K, be fixed Hilbert spaces called the initial and internal spaces, respectively. L 2 (I; JC) = L 2 (I) <£> JC is the Hilbert space 
of square- integrable AA valued functions on an interval / C MA. The (Bose) Fock space over L 2 (I;JC) is denoted by 
r + (/;/C). In the following, a quantum stochastic process will be understood as an operator-valued family (X t )t>o on Hq 



®r + (]R + ,/C). The initial value of the process Xq will be taken as an element xq € B(Ho). The noise space admits the 
following continuous tensor-product decomposition, 



e=0 



r+(R+ /C) =r + ([0,t];/C)®r + ((i,co);/C), (3.1) 

which allows the introduction of a time filtration. In particular, a process (X t ) t>0 is adapted if, for each t > 0, X t (p = ip, 
for all ip <G T + ((t, oo) ; JC). 

Let ^ (/) denote the exponential vector with test function / G L 2 (B. + ,1C) and for some subset S let EXP (<S) be the span 
of exponential vectors with test functions in S; the conserver A (M), creator (g) and destroyer A(g) are operators on 
r + (R + ,/C) defined by their actions on the span of exponential vectors, 

A(M)*(/) : =-i^(e ieM f) 

A (<?)*(/) : =£-*(f + eg) , 
9e e=o 

At(ff)*(/) : = <<?> />*(<?)• (3.2) 
The following notation will be used here 

A] 1 := A (P [0jt] ) , (.g) := # (x [0 , t] ^) , A 01 (g) := A ( X[0 , 4] W) , := t, (3.3) 

for ( > 0, j £ Lj^ c (M + ;/C) and -P[o,t] : / ^ X[o l]®/- They are called the conserver, creator, destroyer and standard 
time process, respectively. One may write Aj 1 , A\° , A 1 } 1 , A® for the creator and destroyer when the intensity g = 1. For 
general g one then has A™ (g) = f* g (s) dA 1 ® and A® 1 (g) = A\ Q (g)\ 

Let I? be a linear domain in T-Lq and let S be a dense linear manifold (called an admissible space) in L 2 (R + ; K) closed 
under the action of the projections -P[o,t], it > 0). An adapted process (X t ) t>0 is said to be based on (T>,S) if each X t , 
is the ampliation to T> ® EXP (Pr 0) t]<SJ <E) T+ ((t, oo) ; /C)of an operator on the domain T> (£) EXP (P[ 0jt ]5) . 
Hudson and Parthasarathy [12] define stochastic integrals of the type 



X t : =x + / {X? ® dA? + X™ ® dA™ + X? 1 ® dA™ + X™ ® dA™ } 
./o 

= a-o+ / Xf®dAf, (3.4) 
Jo 

where the ^A t Q/3 ^ arc adapted quantum stochastic processes based on (T>, S) which are weakly measurable and satisfy 
the local square-integrability conditions below, 

ds\f{s)\ ||A s n u®*(/)|| 2 < oo, f ds\\X^um (f)\\ 2 <oo, (3.5) 

Jo 

for arbitrary t > 0, u <E I? and f € S. 

The notation in (3.4) disguises the fact that the integrators commute with the integrands (when adapted). The algebraic 
manipulation of stochastic integrals is then summarized by the quantum Ito formula for adapted quantum stochastic 
processes X t and Y t , 

d(X t Y t ) = X t dY t + dX t Yt + dX t dY t (3.6) 
where the last term is the quantum Ito correction and can be evaluated using the quantum Ito table 

dA^dAl = dAf , all other products vanishing. (3-7) 

DEFINITION 3.1. Let X t and Y t be adapted processes. One defines Stratonovich quantum stochastic integrals as 

T t dY t := f (x t + \dXA dY t - j dX t Y t := f dX t \ Y t + ^dY t j (3.8) 



The product formula is then X t Y t = xq-i/q + J Q (dX s Y s + X s dY s ^. The integrands X t , Y t do not commute with the 
differentials dX t ,dY t under the integral sign. This lack of commutativity will result in Stratonovich integrals having 
left and right hand forms. To work out general transformations between the two calculi for functions of a process, the 
following lemma will be useful. It is easily established by induction. 

LEMMA 3.2. If X t satisjics the QSDE dX t = X? P ® dA^ then, for polynomial / = f(x) = f n x n , the quantum Ito 
calculus gives 

df (X t ) = f(X t + dX t ) - f (X t ) = f (X t ) af3 ® dAf, (3.9) 



where 



f(X t ) n = f{X t + Xl l )-f{X t ), 

n 

f( Xt ) 10 = Y,fnT,( x ^ + x n pl ^°(x t y 



n Pi,P2 



f (x t r = ^/„^(i t rir(i t +i t T. 

n pi ,P2 
n 

n 

+£/« £ (x t rx^(x t +xiTxi°(x t r, (3.10) 



Pl,P2,P3 

Now, from the definition / (X t )dA^ = [/ (X t ) + |df (X t )\ dA"^ and with the use of the quantum Ito table, one has 

JVQdA? = \[f{X t +X?)+f{X t )]dJ^ + \f{.X t ) 01 dJ%\ 
WQdA? = \[f(X t +X?)+f(X t )]dAl° + ±f(X t ) 01 dA™, 



f{X t )dAf = f(X t )dA 01 



t ! 

f{X t )dAf = f(X t ) 00 dA° t . (3.11) 
Likewise 



dAl'JixT) = i[/(X t + ^ 11 )+/(X t )]d^ t 11 + i/(X f ) 01 dA° 1 , 

dAl°J{x7) = f(X t )dAl°, 

dAfHXT) = \[f{X t + X^)+f(X t )]dAf + l -f{X t f Q dAf, 

dAfJiXl) = f(X t ) 00 dA° t . (3.11) 

Remember though that, as in the classical situation, these notations apply only under the integral sign and are strictly 
nonassociative; that is, X t Y t dZ t ^ X t Y t dZ t . 

THEOREM 3.3. The following transformations exist between QSDEs where the coefficients are related according to 
(2.11) and (2.12), respectively, 

(i) dX t = v (X t ,t)dt + a (X t ,t) {Al° + dA] } 

= d(X t ,t)dt + a(X t ,t){A 1 t ° + dAl } 

= dtv{X t ,t)dt + {Af + dAf } a(X t ,t), 

(ii) dX t = v(X t ,t)dt + fj,(X t ,t) {dA\ l + dA\° + dA\° + dAf } 

= v{X u t)dt + fl (X t ,t) {dAl 1 + dA]° + dAf + dA° t } 

= dt v(X t ,t) + { dA] 1 + dA) + dA]° + dA 00 } fj,(X t ,t) 



Proof. In (i), the process A\° + A® 1 in the Fock vacuum state, corresponds to the Wiener process. Here X\ x = 0, 
X? 1 = X} Q = a (X t ,t) and X?° = v (X t ,t). From (3.11) one has that the reordering of dA\° will lead to the Stratonovich 
QSDE with integrators to the right being related to the Ito QSDE by 

WQjdA™ = f (X t ) dA\° + l -Y. /»" ~° ( X ^) dA t°> ( 3 - 13 ) 

n 

the last term is clearly just \ f (Xt) a (Xt, t) dt and so 

a(X u t)dA\° = ~o (X t ,t) dAl° + (X t ,t) & (X t , t) dA° t ° (3.14) 

This establishes the relationship. The alternative Stratonovich QSDE where the noise coefficients lie to the right of the 
noise is handled with (3.12). 

In (ii), the process A]: 1 + A] + A\° + A® , in the Fock vacuum state, corresponds to the Poisson process. Here X} 1 = 
X}° = X}° = X? — v(X t ,t) = Jj,(Xt,t). To obtain the Stratonovich QSDE with integrators on the right, the terms 
leading to differences are those in dA 11 and dA 10 ; from (3.11) again one has 

J(X7) {dAl 1 + A] } = l -{f(X t + {L) + f(X t )}{dA\ l +A\°} (3.15) 

1 ™ 

n pi,P2 

The coefficient of the last term sums to | {/ (X t + Ji (X t , t)) — / (Xt)}, and so one concludes 

JJXt){dA 1 t 1 + dAf + dA] + dA° t } 

+\ {/ (Xt + n (Xt,t)) f (X t )} {dA\ l + dA}° + dA] + dA™} . (1) 



4 Quantum white noise representation 

There is a more natural way to look at this problem. Take as admissible space the subset S C L 2 (M. + ,JC) got by taking 
the sup- norm completion of the square- integrable step functions on R + . That is, S is the set of square- integrable regulated 
functions and for f € S one is guaranteed that the past and future instant limits f(t~) and f(t + ) exist at each t > 0, cf. 
Dicudonnc [21]. Note that the projection requirement of admissible spaces rules out the Schwartz functions; the space of 
test functions S is in a sense the most natural choice of it as the widest space on which integral approximations can be 
based and also contains the functions of bounded variation which are the natural space on which to discuss functional 
integral transforms. One defines functionals S± (t) and 5* (t) on iS by 

(S± (t) ,f)=f (t±) , (S* (t) , /) = \{f (t+) + f (t- )}. (4.1) 

That is, the action of (5* (t) on a function of a time variable is to give the average of the immediate past and future values 
at time t. The action of these functionals can be extended by linearity from S to S' if the following identifications are 
made: 

(M*)>M*)> : =M*-«), 

(S+(t),6.(s)) : =S-(t-a), 
(S + (t),S+(s)) : =5*(t-s), 

(S.(t),S.(s)) : =S.(t-a), (4.2) 
DEFINITION 4.1. A basic step function on R™ is the characteristic function of a set of the form {(ti, ■ ■ ■ ,t n ) : < tjm — c\ < ■ ■ ■ < 
where 1 < r < n, the j(l), ...,j(r) are distinct elements of {1, n}, and cq,Ci, ■ ■ ■ ,c r are constants, and also where any 



of the strict inequalities can be replaced by ordinary ones. A simple function on R™ is a finite linear combination of such 
step functions. The space of multi-dimensional regulated functions iS„, is taken to be the sup- norm completion of these 
simple functions. 

For / G S, the individual limits of / (t\, ■ ■ ■ ,t n ) exist as the tj — > a^ under the conditions t a (i) < 



• < t a (n) for a11 



a,- G 



and for each permutation a of {1, 



,n}. 



DEFINITION 4.2. Quantum white noises are defined on EXP(S) by J± (t) 

Explicitly for / G <S, one has a± (t) *(/) = / (i±) * (/) ; e tc. 

From (4.2), the nontrivial commutations between the a\ (t) are given by 

[a-(t),al(s)} = 6+(t-s), 

[«+(*) ><*-(«)] = 5-(t-s), 
[a+ (t) , al («)] = 5* (t - s) , 
[a_(i),a f _(s)] = 5*(t-s), 



A* (S ± (t)) and a! (t) := A» (5* (i)). 



(4.3) 



The linearity of the extended functionals then implies that the pair of processes |a| (t);t > Oj satisfy the following 
canonical commutation relations (CCR) 



[a, (t) , a, («)] = = [at ( t ) , a t ( s )], [ a>> ( t ) , a t ( s )] = s ) 
One further has the following functional distribution 

8°# {/(*)<*!(*) + /(*)* «!(*)} 



cxp 



exp 



--II/II 2 
2 IU 11 



(4.4) 



(4.5) 



where / G 6>, and ||/|| 2 = / °° \ f (t)\ 2 dt, and the expectation is in the Fock vacuum state VP (0). The key feature of (4.4) 
is that, along with the specification of the state (4.5), it contains all information concerning the chaotic expansions. Any 
integral of the form J" R „ dt\ ■ ■ ■ dt n ip (ii, • • • , t n ) at^ (ii) • • • at^ (t n ), for ip G S, can be evaluated and, in particular, one 
may take p to be simplicial. 

The following connection [10,17] exists between the quantum stochastic calculus and white noise calculus. 
THEOREM 4.3. Let (X t ) t >o be the solution of the QSDE dX t = X" <E> dA" with the (x^ ^ adapted processes 
based on (D,S), then the QSDE can be represented as 

dX t = {al(t)Xl 1 a.(t)+at(t)Xt° + X° 1 a m (t)+X? > }dt 

= M(*)] 



xl° 



x? 1 
x} 1 



1 

a* (t) 



Remarks 

i) The Wiener and Poisson processes are represented by B t = J * (at (s) + a* (s)J ds, N t = / * (1 + a* (s))* (1 + a* (s)) ds, 
respectively. Their chaotic expansions can readily be obtained. 

ii) A Stratonovich QSDE is an equation of the form 



dX t 
dX t 



{4 {£) a* (t) E t R + at (t) F t fl + a, (t) G t * + H*} dt 
\W t L al (t) a, (t) + T t L at (t) + ~G t L a* (t) +~H~ t L } dt 



(4.7a) 
(4.7b) 



Equation (4.7a) is the left handed version, and (4.7b) is the right handed version. 

If the conserver terms are ignored, then one sees that the quantum Ito calculus corresponds to the Wick ordering scheme 
while the quantum Stratonovich corresponds to the Wcyl scheme. This point of view can be of help in understanding 
the related distinctions which arise in the theory of phase space path integrals [22]. The anti-Wick ordering scheme gives 
time-reversed quantum Brownian motion, 
iii) The classical Wiener integral is then represented as 



X t = X a + J {v (X„s) + [at (s) + a, (s)] o (X s ,s)}ds 



(4.8a) 



X t =X Q + f (X s , a) + aj (a) a (X s , s) + a (X s , s) g, (a)} 



da 



where v,a and v,d are related by (2.11). 

Likewise the classical Poissonian integral is represented as 



X t =X + f {v (X s , s) + [l + a, (s)] f [1 + a, («)] /x (X„ «)} ds (2) 



or 



X t = X + J { v (X„ s) + [1 + a. (s)] f A (X s , s) [1 + a. (a)] } , 
where v,/x and £>,/2 are related by (2.12). 

iv) Let X t = xq + J Q X"* 3 <g> dA" 13 ,Y t = y + f Q Y" 13 ® gL4"^ be processes based on (2?, <S) , then 

XldK, = J^ds(^x a + J ds {at (s) X} l a, (*) + oj (a) X s 10 + X s 01 a* (s) + X s 00 } 
x {4 (a) n 1 V (a) + at (a) Y s w + Y s 01 a* (a) + Y s 00 } 

ds {at (s) X s Y s n a* (s) + a\ (s) X S Y S 10 + X s Y s 01 a, (s) + X S Y S 00 } 



+ J ds J du5* (u - s) [at (tt) X^Y^a* (s) 

+at(u) X?Y, W + X°u% n a* 00 + KV] 

= f X s dY s + \f ds [Xf + at (s) X s u ] [F s 10 + Y s n a* (a)] . 
Jo z Jo 

That is, Jl X~ s dY a = J* X s dY s + \ f* dX s dY s . 

Thus the algebraic product X t Y t of the white noise representations gives the correct product as quantum stochastic 
processes. 

v) From the above considerations, it follows that the following formal manipulations are allowed 

T r T 



f X t ®dAf-( dAfX t = f [X t ,a*{t)]dt 
Jo Jo Jo 



dt / ds [at (t) X^a, (i) + at (t) X™ + X^a, (t) + X?°, a, (a)] 
o Jo 



[ dt [ dsS, (t - s) [X^a, (a) + X} ] 
Jo Jo 

dt [X^a, (t) + Xl°] 
dt [X] 1 <g> dAj° + X™ ® dA° t 



1 " T 



2 
1 

2 Jo 



1 " T 



More generally, for polynomial /, 

df (X t ) = {at (t) f(X t ) n a* (t) + 4 (t) f(X t ) w + f(X t ) 01 a* (t) + f(X t ) 00 } dt, 
and formal manipulation leads to 

[/ (X t ) , a, {t)] = -\f{X t )^a» (t) - \f(X t ) w , 

or under rearrangement 

l /(X t +X t n )+/(X t ) q. (t) = a, (t) J{X7) - \f (X t ) 10 . 
This is the same as the third relation of (3.12). 



5 Asymptotic quantum stochastic limits; convergence ansatz 

Let /Co be a fixed Hilbcrt space and S(t) = e tM be a strongly continuous one-parameter unitary group on ICq. If JCq is 
the subspace of JCq such that dt \ {S (t) g, f)\ < oo whenever /, g G /Co, then introducing the sesquilinear form 

/>oo 

7 (/,<?)=/ dt(S(t)g,f) (5.1) 



the Hilbcrt space completion of /Co, with subspace {k : 7(fc, fc) = 0} factored out, will be denoted by /C. The inner 
product on /C will be taken to be that inherited from /Co which will be denoted (•, -)^. For each g G /C the mapping 
t S (t) g is Bochner intcgrablc. The following operator is well defined on JC: 

S + :=l° S(t)dt^^ = 7 r5(n) + l Pvl (5.2) 

and the notation S- := T = 2i?eS+ and £ := ImS+ shall be adopted. Noting that T > and that [T, E] = 0, it 

follows that there exists Z = X + iY with X and Y self-adjoint on JC such that 

X 2 = -T, XY + YX = -S. (5.3) 

For each 7 > 0, one considers Ha a copy of T + (JC) with Fock vacuum vector denoted as "J a and A^ x (•) the creator/destroyer 
maps. For g G /C,t > 0, one introduces the following operators on "Ha 



o* (t,ff) = A« Qs(t/A 2 ). ) 



(5.4) 



From the CCR one has 

oa (t,g) ,a[ (sj) = -^(S f-^F") 9>f)z ( 5 - 5 ) 

which says roughly that these operators have auto-correlation time of the order A 2 . As A — > 0, one expects these operators 
to become white noises, however, taking account of the previous section, the functional limit of (5.5) can be interpreted 
more concisely as 

[a (t, g) , 0) ( s , /)] = (S+g, f) K S+ (t - s) + (S.g, f) K S- (t - s) (5.6) 
The limit operators a" (t,g) are interpreted as the operators defined on H := T + (L 2 (R)®/C) by 

a (t, g) * (0®/) := (6+ (t) ®Zg + <5_ (t) ®Z*g, fef) <J (0®/) (5.7) 

with a'(t,g) the adjoint of a(t,g). 

THEOREM 5.1. The processes {a\(t,g) : t > 0,g G JC} on 1~L\ with the state \Pa converge in Fock vacuum expectation 
as A — > to the quantum white noises {a?(t, g) : t > 0,.g G JC} on H with the state \f f . That is, for each integer n > 0, 
the following limits hold on S n , 

ton<* A ,af > (ti, 5l ) • • -af l) (*„,<?„) v& A )« A = (^,a^ (t u9l ) ■ ■ ■ a«"> (t n ,g n ) (5.8) 
for all ti, • • • , t„ > 0, t/i, (/n S /C and choices of creators and/or destroyers. 

Proof. From the Gaussianity of the pre-limit processes it suffices to consider the two-point functions. For <j) g S2, let 

h = I dsdt c/}(t,s) -j(S f^T2~) g,f)tc 

JR+XR+ A V A / 

00 />u/A 2 

du / dr <f> (u + X 2 t,u - A 2 t) (S (t) g, f) K , 
J-u/X 2 

where the change of variables u = t + s, r = (t—s) /A 2 was made. If Jg is the expression obtained from I\, by replacing 
the i-limits of integration by ±00, then \I\ — J\\ — > as A — > + uniformly since /, g G JC. Moreover, since (f> is 
L 2 -regulatcd, J\ converges uniformly to 



du {(f>(u + ,u )(S+g,f)ic + ^(u ,u + ) (S-g,f)ic} 



which establishes (5.6). 

Next the CCR of the limit noises must be checked: 



[a(t,g),aHsJ)] = (8+ (t) ®Zg + S. (t) ®Z^ g,S+ (s) ®Zf + 5- (s) ®tf f) n 
= ((Z^ 2 + ^Z + ZZ^)g,f) K 5 + (t-s) 

+ ((Z 2 + h&Z + Ztf)g, f) K S- (t-s), 
but (Z+) 2 + \{&Z + Ztf) = 2A 2 - i{XY + YX) = \T + zS = S+ and so (5.6) is recovered. ■ 

DEFINITION 5.2. Let (h) denote the exponential vector on H\ = F + (/C) for k e JC. Collective exponential vectors 
are defined, for / e S, by 

(5.9) 



* (/, k, A) := * A U °° dt f (t) S (i/A 2 ) fc 



The set of such collective exponential vectors is denoted EXP\(S, JC). 

Noting that (\& (/, k, A) , ^ (/', k', — > exp{7(fc,fc') (/, /')l 2 (k+) } as A — > + , it is natural to associate the limit 
* (f®r?k \ in U with * (/, k, A) by virtue that 7 (fc, fc) = (rfc, k) K = ||rV2fe||^.. 

In the following we shall understand all processes to be over a common domain D of Ho, as outlined in Section 3. A 
family of operators (Z t (A)) on Ti.Q§fH\ is said to converge to a process Z t on 7io§[H as A — > + weakly in matrix 
elements [9] if for all n, ti, ■ ■ ■ ,t n , and for all u, u' £ T>, f,f € S and k,k' € K, one has 

lim (u ® * (/, fc, A) , Z ti (A) ■ • ■ Z tn (A) u' <g> * (/', k', A)) = (« ® * , Z ti (A) ■ ■ • Z tn (A) «' ® * (/W)) . 

Typically one would also like Z t (A) to converge to an adapted process and for this reason the next definition is formulated. 
DEFINITION 5.3. A process (Z t (A)) t>0 on W §M is said to be adaptable if, for each t > and for all s > t, g e JC, 
one has 

Z t (X),a{(s,g)\=O ttS (X), (5.10) 

whereby a process W = Wt. s {h) is said to be Ot >s (A) if J Q T ds J Q S dtWt, s (A)vanishes in matrix element limits. Further one 
requires that, for all u e Ho, f S <S and g G JC, and A in a neighborhood of + , 

t 

ds \\Z S (A) u ® * (/, fc, A)|| 2 < C, (5.11) 

Jo 

where C depends at most on u, f and g. 

Let A t Q/3 (A) be operators on Ho®H which depend on a\ (t, •) for s < t. Note that, since the a\ (s,g) have finite auto- 
correlation time, it does not follow that they will commute with X"^ (A) for s > t, however the assumption shall be 
made that they are adaptable. Consider then ODES of the type 



d 
It 

Xo (A) = x £B (Ho) 



-X t (X) = a[(t,g 1 )Xl 1 (\)a x (t,g 2 ) + a[(t,g 3 )Xl°(\)+X° 1 (X)a x (t,g 4 )+X^(X), (5.12) 



In this form the ODE has creators and destroyers in normal ordered from. The objective of the remainder of this 
section is to show that (5.12) converges to a well-defined QSDE under the ansatz that its coefficients are adaptable. For 
convenience, the K-state will be a fixed g and it will be supposed that (S±g,g)ic = \ for the rest of this section. The 
g-dependence in the prc-limit and limit noises and collective exponential vectors will be dropped, and one notes that in 
this case the limit noise is just the fields a* (t) introduced in the previous section. 

THEOREM 5.4. Suppose At (A) is the unique solution to a Wick ordered ODE with adaptable coeficients A t (A)) 
converging to processes A"' 3 so that the QSDE dXt = A t a/3 ® dAf , Xo = xo, has again unique solution At. Then At (A) 
converges to At , weakly in matrix elements as A — > + . 



Proof. Part (i). Note that one has now taken j(g,g) = 1 and writes ^ (/) for $ (ftgiTzg^. Let (pi,(p 2 € ~Hq, 
f,ke L 2 (M+), then as A -> 0, 



fo®* (/, A) , {X t (A) - x } V2 m> (k, A)) 

i 



o 



o 



A 



A J 



fc(v) 



g)X 1 s 1 (X) I dv-^(g,S 



A 2 



9) 



d S <^, {/ (*)* X^k (s) + f (S)* X™ + Xfk (s) + X° } V2 ) Ho 
(f) , f Xf ® dAf <p 2 m> (k)}u ®n ■ 



This amounts to saying that the approximation holds at the level of the QSDE, that is, Xt (A) converges to Xt in first 
moment. 

Part (ii). It is easy to see that higher moments X tl (A) , • ■ • ,X tn (A) decouple in collective coherent states, under the 
assumptions of the theorem, provided the tj are distinct. To deal with the case where several indices are equal, one 
actually shows that such moments converge to the appropriate equal time Ito products. This is established by first 
showing that if Y t (h) is a similarly described process, then the product Y t (h) X t (h) converges to X t Y t in QSDE. 
Multiplying the pre-limit operators, one obtains 

(X t (A) - x ) (Y t (A) - y ) 



du 



dv {a[ (u) Xl 1 (A) a x (u) + ■ ■ ■ } {a[ (v) YX] 1 (A) a x («) + ••■} , 



there are 4x4 separate integrals here, examination of the first will be sufficient to determine the general pattern; one 
has 



du / dv 4 («) K 1 (A) a x (u) a[ (v) Y v n (A) a x (v) 
Jo 

f duj^ dval(u)Xl 1 (X)^al(v)a x (u) + ^(g,S^^^g)^Y v 11 (X)a x (v). 



If v < u, then a x (u) can be commuted with Y^ 1 (A) up to an Uj „ (A) error. Similarly, if u < v, then 4 (v) can be 
commuted with X^ 1 (A) up to an O u>v (A) error, thus the above equals 



du 4 (u) Xl 1 (A) | J dva{ (v) Y v n (A) a x (u)j a x (v) 
dua\ (u)X^ (A) |y dvO v>u (X)a x a x (v) 



+ dval (v) \ / dual (u) Xl 1 (A) a x (u) \ Y v n (A) a x (v) 



dv / du 4 (v) a x (u) Y^ 1 (A) a x (u) O u<v (A) 



du j dv (g, S 



A 2 



In the limit A — > leads to 











/ n" 


<g> dA" 


/o 


JO 








f-U 


+■ r ^ 




/ n" 


^0 




/o 



ff )4 (A)y^(A)a A («) 



dA 1 



The last term is an Ito correction. Such a term results whenever the pre-limit term 4 ( v ) a x ( u ) is present and put into 
normal order. There are 4 such terms and they lead to the usual quantum Ito correction 



lim (X t (A) - x ) (Y t (A) - y ) = (X t - x ) (Y t - y ) 

A->0+ 



Thus the Ito calculus is picked up in the limit and the convergence for all moments can be derived by induction. ■ 
Here we have not attempted a most general statement which might be formulated by defined processes as equivalence 
classes of sesquilinear forms on the appropriate space of exponential vectors. Instead, we assumed existence and unique- 
ness of solutions to begin with, however, situations where this can be established will be presented in the next section. 
THEOREM 5.5. Let X t (A) = x + du {a{ (u) X^ 1 (A) a x (u)-\ j with adaptable coefficients (A) such that 

X t (A) , X" (A) converge respectively in weak matrix elements to X t ,X" (A) (uniformly for mixed matrix elements) us 
above, then the following limits hold us quantum Stratonovich integrals: 



ds a\ (s)a x (s)X s (A) 



dsX s (A) a\(s) a x (s) 



ds a\ (s) X s (A) 



ds X s (A) a\ (s) 



dsX s (A) a\ (s) 
ds a{ (s)X s (A) 



dA x }X 



X»dAi 



dA°}X. 



X,dA 



10 



01 



X H dA 



dAl°X 



(5.13) 



Proof. The treatments of these limits are very similar, each involves at most one reordering, and it is enough to work 
through just one of them: 

dsa x (s)X s (X) = [ dsX s (\)a x {s) 



+ ds du [ax (s) , a[ (u) X^ 1 (A) ax (u) + ■■■'_ 



j ds X s (A) ax (s) + j ds j du—(g,S 
Jo Jo Jo A 

+ / ds du a\ (u) Oxax (u) . 



g){Xl 1 (X)ax H+^°(A)} 



Now, for Rx(u) the matrix elements of an adaptable process between arbitrary collective exponential states, one has 



lim ds du —7r(q,S 
a^o+ J J A 2 U ' 



g)Rx (u) 



lim 

A-J-0+ 



ds 



s/X 2 



dr (g,S(r)g)Rx (s-r/X 2 ) 



R (s) ds, 



where R(u) is the associated limit; from the standard notion of weak convergence in matrix elements from T-Lq^Hx to 
Ho^T-L one sees that 



lim 

A-S-0+ 



dsax (s)X s (X) 



X, <g> dA 



dA^X s . 



01 



1 

2 7 



li 



<dA 



10 



X 



10 



dA 



O0\ 



6 Asymptotic quantum stochastic limits. Uniformly convergent situa- 
tions 

The simplest way to bypass the adaptability ansatz introduced in the previous section is to consider only linear ODES. 
It is possible, however, to construct genuinely nonlinear examples starting from the linear case. This, in fact, is very 



natural: all important dynamical equations in science are nonlinear with the sole exception of the Schrodinger equation, 
this however leads to the Heisenberg evolution which is generally nonlinear. The program of this section is as follows. It 
is first of all shown that solutions U t (A) to linear Wick-ordered ODES converge without adopting an adaptability ansatz. 
If the limit process Ut, is unitary, then it is shown that processes X t = u}x$Ut, for some xq G B (Hq), satisfy nonlinear 
QSDEs. In a sense, this class of solutions is the most interesting as they provide examples of quantum dynamical variables 
having a symplectic evolution. 

THEOREM 6.1. Let x , C a p G B (Ho) for a,/3 G {0, 1} and let X t (A) be the solution to the linear ODE 

X t (A) = a{ (t) C u X t (A) a x (t) + a{ (t) C w X t (A) + C m X t (A) a x (t) + C m X t (A) (6.1) 
with Xq (A) = xq, then for all H.§ and f,h £ S 

lim (u®* (/, A) , X t (A) (fc, \))u ®-H = (um (/) , X t vm> [k))n ®H, 

X— >0+ — — 



where X t is the solution to the QSDE 



dX t = J2 C ^ Xt ® dA t P > X ° = x 0- ( 6 - 2 ) 

a,/3 



Proof. Let a" (t) = a\ (t) for a = 1, and = 1 for a = 0. Eq. (6.1) can be re-written as 

Xt(X) = Y,a a x (t^C aP X t (X)a f i (t) . 

a,/3 

As this is linear, the formal iterative series expansion exists, 

oo 

X t (A) = x + Y, x t n) (A) 

n=l 

where 

A t (n) (A) = E / dti-'-d^ 

«i^i,-,«»,?„ i>t " > "- >il - 

x (a^ 1 (ti) • • • a°" (t„)) f C aifll • • • C an/ j B a:o a^ 1 (*i) ' • ' a i" (*n) 

In particular, 



A),X 4 (A) «2®*(fc,A))-H ®H = V] / dti---dt n {u,C ai p 1 ---C an p n XQv) Ho 

oi,fl I .-.««»./9n" / * >t " > '" >tl - 

X/A 1 (*l)*-/r , (*n)**A 1 (*l)-*A B (*») ( 6 - 3 ) 

where /" (t) := /J" f (t + \ 2 u) (S (u) g, g)du for a = 1, and = 1 for a = 0, etc. 

Expression (6.3) is bounded by Co ||u|| |H| (c(Tg, g) WfW^ \\k\\ 00 t) n where c := max |C Q( gj|. The scries expansion is there- 
fore uniformly convergent, with (6.3) converging to 



y~] / dt 1 ---dt n {u,C ai p 1 ---C an p n x v)u 



X /«i (*!)*..•/«» (t„)*AA (ti)...^- (<„), 

where J 1 (t) := (S+g,g)f (t)"and /° (i) := 1. Resumming gives the correct matrix element for the process X< described 
in the statement of the theorem. ■ 

The expediency of Theorem 6.1 compared to those employed in the papers of Accardi, Frigerio and Lu [9] comes about 
from the fact that here the limit is anticipated via white noise operators. The mechanism is transparent because the 
pre-limit and limit representations are Wick ordered. The difficulty encountered there was that the prelimit fields were 
Weyl ordered (which is generally the ordering natural to equations of elementary physics) and so enormous efforts were 



spent in re-ordering and the subsequent identification and treatment of negligible terms. Once it is known that the 
process converges, however, the adaptability ansatz can be dispensed with for the coefficients of the ODE. Part (ii) of 
Theorem 5.4 affirms that the multi-moment convergence of the process occurs. Therefore the following conclusion is 
reached. 

THEOREM 6.2. Let X t (A) be the solution of the linear ODE (6.1) then X t (A) converges weakly to the solution of 
the QSDE (6.2) in matrix elements. 

The problem of obtaining a unitary process Ut obeying a linear QSDE has been tackled [12, 171 and can be summarized 
in the next theorem. Here a" (t) = as(t,g) with (S+g,g) taken as n = 57 + icr. The CCR for the noise is then 
[a(t),a)(s)] = k5 + (t-s) + kJ_ (t - s). With the identifications dA] 1 = a) (t) a (t) dt, dA]° = a) (t) dt, dAf = a it) dt 
and dAf ° = dt one is led to the nontrivial component of the lto table being dA^dA] = 7cL4"' 3 . 

THEOREM 6.3. The general unitary process [23] (U t ) t>Q driven by white noise processes a" (t) is given by (Uo = 1) 

1 



U 



1, 



n 



« T (*) 



-iH-WL -tfW 



L 



W-l 







u t 





-i [Eat (t) a (t) + FcS (t) + F^a (t) + G] 



(6.4) 
(6.5) 



where E, H, G are self-adjoint and W is unitary on Ho, and the coeficients of the lto QSDE (6.4) are related to those of 
the Stratonovich QSDE (6.5) by the relations 



W 



l-in*E 
1 + ikE ' 



(1 + inE)~ l F, H = G + F^ 



a - \nf E 
1 + |k| 2 ^ 



-F. 



(6.6) 



One notes that the Stratonovich QSDE takes the form U t = -iT t U t , where T t = Ea) (t) a (t) + Fa) (t) + F^a (t) + G 
can be interpreted as a quantum stochastic Hamiltonian. 

THEOREM 6.4. Let E,F,G e B(H ) with E and G self-adjoint then the family of processes U\ (t) satisfying the 
ODES 

Ut (A) = —iTt (A) U t (A) , U (X) = 1, 

with 

T t (A) = Ea\ (A) a t (A) + Fa\ (A) + F^a t (A) + G, 

converges weakly in matrix elements to the unitary quantum stochastic process described in Theorem 6.3. 
Proof. One begins by noting that 



[a x (*) , U t (A)] = -1 



a x (t), / dsT s (X)U s (X) 



[* 1 
Jo ^ 



t-3 



g) {Ea s (A) + F} U s (A) 



= -i(S+g,g) {Ea t (A) + F} U t (A) + O (A) . 
Rearranging the ODE from Wcyl to Wick ordered form gives 

Ut (A) = -a[(t)(W-l)U t a x (t) + ^-a{(t)LUt 



(6.7) 



7 



n 



V7 



tfWUtax (t) - f -Dh + iH)U t + (A) 



one recognizes the prc-limit form of (6.4) and the results of the previous section imply convergence. 



The weak convergence in matrix elements, however, gives much more. 

THEOREM 6.5. Let x a € B{H ), and set X t (A) = \j\ (X) x U t (A), then X t (A) converges weakly in matrix elements 
to the quantum stochastic process X t = UjxoU t - 



Given X t — lllxoUt, then one has the stochastic evolution equations 



X t = (l}x U t + Ujx U t 



1 



[x t ,e t 



(6.8) 



where Q t '■= UjTUf. The second part takes the form of a stochastic Heisenberg equation. It is a straightforward 
calculation to show, either by converting to the Hudson-Parthasarathy calculus using the quantum Ito formula and 
inverting back or by the now standard manipulations using the white noise calculus to put to Wick order, that the QSDE 
for X t , is 



X t = 



-a) (t) (w}X t W t -X^a(t) 



:<zt (t)W}[X u L t ] 



X t ,L t 



W t U t ax (*) 



L f , Xt 



L t --L\ [Xt,Lt]-i[X tt H t ] 



(6.9) 



The Stratonovich version of this will now be computed. 

LEMMA 6.6. Let Ut be the solution to the QSDEs (6.4) and (6.5), then one can perform the following commutations 
under the integral: 



uW (*) 
Uja (t) 



{a 1 (t) u\ + ik*U^F^\ (1 - iK'E)' 1 , 
la (t) u\ - IkUIf} (1 + inE)' 1 , 

U}a){t)a(t) = {a? (t) a (t) Xj\ — iko) it) \j\ (1 — in*E)~ l F 

+in*a (t) Uj (1 + inE)" 1 F f + \h\ 2 U}f (1 - in* Ey 1 F^ 
+\k\ 2 u\f^ (1 - i^Ey 1 F} (1 - 2oEy X . 

Proof. Before demonstrating the proof it is important to point out that (6.12) cannot be obtained from (6.10) and (6.11); 
this is due to the afore-mentioned nonassociativity (that is, the product of Stratonovich integrands is not defined as the 
Stratonovich integrand of the products). To be consistent, the left-hand side of (6.12) should be written as U/at(t)a(t); the 
symbol U/at(t)Z(t) has not been defined. However, the underbar and overbar notation has been dropped for convenience. 
To prove (6.10) one has directly 



[a(t),U t ] = 



[a (t) , (Ea) (t) a (t) + Fa) (t) + F*a (t) + G) U s ] ds 



(Ea (s) + F) {kS+ (t - s) + kJ- (t - s)} U s 
= -iK(Ea{s) + F)U t . 
The right-hand side is not properly normal ordered; however, it can be rearranged to give 

a (t) U t = (l + inEy 1 {U t a (t) - inFU t } . 



.13) 



Note that this can be deduced from the pre-limit expressions and corresponds then to (6.7). Eq. (6.10) is the conjugate 
equation to (6.13). The derivation of (6.11) is similar. 
To derive (6.12), note first of all that 

[a f (t) a (t) , T t ] = Ea) (t) a (s) {kS+ (t-s) + K*S- (t - s)} - Ea) (s) a (t) {nS+ (s -t) + n*5- {s - t)} 

+Fa) (t) {kS + (t - s) + kJ- (t - s)} - F^a (t) {k5+ (s-t) + n.J_ (s - t)} . (6.14) 



Therefore, 



[a f (t) a (t) , U t ] 



[a f (t)a(t),T s U s ] ds 
-iEa^ (t) a (t) (« - «*) U t - iKFa) (t) U t + in* F^ a (t) U t . 



(6.15) 



Again (6.15) is not wholly ordered, however, whereas the manipulations are not associative, they are linear and so a" (t) Ut 
can be replaced using (6.13), etc. This will lead to the conjugate equation to (6.12). ■ 
The first form of the Heisenberg equation (6.8) can be written as 

X t = U\ (at ( t ) a (t) - [x , E] + at (t) - [x ,F] + - [x ,F^] a(t) + - [x , G}\ U t . (6.16) 

The conversion to Stratonovich form (left-handed) can now be made using the lemma. One obtains 

X t = a) {t)a{t){l-2aE t )- 1 ± [X t ,E t ] 

+at (t) {-in (1 - iit'Et)- 1 F t (1 - aEt)' 1 - [X t , E t ] + (1 - iK*E t )~ l - [X t , F t ] 

+a (t) !^ik* (1 + inEtY 1 F\ (1 - oEtY 1 - { [X t ,E t ] + (1 + inEt)' 1 - { [x t , F t f ] J 

+{- [X t ,H t ] + \ K \ 2 F t (1 + inEtV 1 F} (1 - 2aE t y 1 - [X t , E t ] 
% i 

+ |k| 2 F\ (1 - in*E t y X F t (1 - 2 ( rS t ) _1 - [X u E t ] 

i 

+iK*F] (1 - i^EtT 1 - [X t ,F t ] - nF t (1 + iKE t y X - \x u F}] . (6.17) 
Eqs. (6.9) and (6.17) are equivalent and give the most general Heisenberg equation for this situation. 
Remarks 

The dynamical evolutions considered in this section are broad enough to include the classical SDEs considered in Sections 
2 and 3. 

Let xo,po be canonically conjugate variables on W . Here the specification k = |, W = 1, L = —ia(xo)po, and 
G = \ (v (x )po +pv {x )) leads to (6.9) being the diffusion described in (4.8b). The QSDE has unbounded operator 
coefficients, but this can be dealt with using standard techniques [24]. Here E = 0,F = p a {x ), and H — G, and so 
(6.17) reduces to 

Here Pt = Ujpollt and since the evolution is unitary, [X t ,Pt] = i for each t > 0. The last term in (6.18) is therefore 
|i [a (X t ) ,Pt] a (Xt) = {X t ) a (X t ) which is the correct drift term, and so (6.17) reduces to (4.8a). The classical 

relations are therefore recovered. 

Likewise, to obtain Eqs. (4.9) for the Poisson driven processes, one sets k = h and takes E = F = F^ = H — 
\{Pqv{xq) + v(x )p ) = ^(poM( x o) + ^( x o)Po)- In this case W can be represented as the change-of-variable operator 

(Wf)(x) == Jfef (u (x)) where u = u(x) := x + p,(x). In this case H^xoVK = x$ + fi (xq), with xq interpreted as the 

multiplication by x operator on T-Lq := 1? (R). 
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